Winding up superfluid in a torus via Bose Einstein condensation 
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We simulate Bose-Einstein condensation (BEC) in a ring employing stochastic Gross-Pitaevskii 
equation and show that cooling through the critical temperature can generate spontaneous quantized 
circulation around the ring of the newborn condensate. Dispersion of the resulting winding numbers 
follows scaling law predicted by the Kibble-Zurek mechanism (KZM). Density growth also exhibits 
scaling behavior consistent with KZM. This paves a way towards experimental verification of KZM 
scalings, and experimental determination of critical exponents for the BEC transition. 



Manifestations of symmetry breaking associated with 
BEC formation are spectacular and diverse pj-[7| . Forma- 
tion of BEC breaks the U{1) symmetry of the phase of 
the condensate wave-function. We show that when BEC 
forms in a ring cooled through the critical temperature, 
this symmetry breaking may result in spontaneous ro- 
tation in the newborn condensate, with long persistence 
due to topological stability of quantized winding number 
W. We demonstrate that resulting winding numbers sta- 
bilize soon after the transition, and the variance of their 
distribution follows the scaling predicted by the Kibble- 
Zurek mechanism (KZM) ^8|, i9|]. Density growth of BEC 
also exhibits behavior consistent with KZM. 

Creation of non-adiabatic excitations (often in form of 
topological defects) when the system is driven through 
the second order phase transition is the focus of KZM. 
In its usual form, KZM predicts density of defects scales 
with the quench rate, with the exponent given by the 
function the critical exponents of the underlying equilib- 
rium phase transition. This has been tested numerically 
in number of studies, see, e.g., [isl, [l6| . However, on 
the experimental side, while the studies to date confirm 
key qualitative predictions of KZM (creation of topologi- 
cal excitations) [l|, Q, [HI , its key quantitative prediction 
(scaling of their density with rate of quench) has not 
yet been convincingly demonstrated (see, however, [l5| 
for suggestive indirect evidence). The difficulty involves 
controlling sufficient range of quench timescales as well 
as counting defects. We study how scaling laws involving 
experimentally accessible quantities, namely, variance of 
W and non-adiabatic response time for density growth, 
pave a way around this long-standing hindrance. 

Recently, persistent circulation of BEC in toroidal trap 
has been achieved experimentally by stirring the cloud 
Q. In particular, with the advent of circular trapping 
potential for BEC [10], the possible experimental test- 
ing of this theory is around the corner. Also, density 
growth in BEC formation with variable cooling rate has 
been already studied [17]. We show that a similar set up 
will allow both verifying KZM for density growth, and 
measuring critical properties of superfluid circulation |8| . 

We consider a BEC in a quasi- ID ring of circumfer- 
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FIG. 1: Sequences of isodensity surfaces for growing condensate 
for a single quench realization with tq = 0.01 for flat (top line, 
time flowing from a to c) and tilted (bottom line, from d to f; tilting 
angle a = 3°) rings. Color represents the phase of the condensate 
along the ring. Formation of random equi-phased patches at the 
onset of condensation is clearly visible. The color shows that the 
magnitude of the winding number W for the condensates in the flat 
(c) and the tilted (f) rings are 1 and 2 respectively. 



ence C, an idealization of quasi- ID toroidal geometry. 



e.g., see 



1 it using the stochastic Gross- 
Pitaevskii equation (SGPE) [if 
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where <p = |^(x)|e*^^^^ is the condensate wavefunction 
and r]{x^ t) is the thermal noise satisfying the fluctuation- 
dissipation relation {r]{x^t)r]^ {x' ^t')) = 2'yT5{x — x')S{t — 
t')^ with 7 representing the dissipation, T the noise tem- 
perature, g the non-linearity parameter and — e the chem- 
ical potential [18]. Leaving aside the noise and dissi- 
pation, the above system can be described by the en- 
ergy functional S = + U{\(l)\)]dx, where 
U{\(j)\) = e|0p + ^^1^1^- Extremizing the energy func- 
tional we obtain (/) = for e > and = ^/\e\/g exp (iO) 
for e < 0, where 6 is the wave function phase (e = is 
the critical point). We induce the phase transition by 
quenching e: 
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from an initial e > to a final e < 0, and allow the 
system enough time to thermalize initially and stabilize 
eventually [20|. The critical point is crossed at = 0- 

When BEC is formed via cooling through the criti- 
cal point, non-adiabaticity is enforced by diverging re- 
laxation time r and healing length ^ Q : 



(3) 



Here v and z are the critical exponents, and tq are 
determined by the microscopic details of the system. Ac- 
cording to KZM, as e approaches 0, after the instant —i 
the relaxation time r(— £) exceeds the timescale |e/e|^^_£ 
of the change imposed by quench, and the state of the 
system freezes. Its order parameter behaves impulsively 
(i.e., remains frozen) within an interval between ±t and 
starts dynamical evolution again thereafter. KZM gives 

r(e(-£)) = |e/e|,=_,-. (4) 

For the linear quench, Eq. (j2j), we get from Eqs. (j3j), (j4]) 
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where ^ = As e becomes negative, condensate starts 
forming with a phase profile 0{x) consisting of random 
patches: phase is approximately uniform over the length- 
scale ^. On each such patch, its value is chosen indepen- 
dently (different stages of condensate formation in a uni- 
form ring are shown in Fig. [1] upper row). The phase of 
each patch is randomly chosen. Therefore, we can esti- 
mate the variance of the total phase Oc = /^^^^ dO{x)^ by 

considering the sum of ~ C/^ uniform random variables, 
each having the same variance 7r^/3 {9 taking any value 
between ±7r). This implies Gaussian distribution for Oc 
(in qualitative agreement with inset of Fig. [3|) with vari- 
ance cr2(6>c) = (Oc) = {7T^/3){C/i). As the wave function 
is single- valued, we must have 6c = 27rVK, where W is the 
integer winding number. So, Eq. (j5j) predicts dispersion: 



ct{W) = ^/W^ 
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using the mean field z = 2 and v = 1/2 [3]. Spatial 
gradient of 0{x) gives local flow velocity. Therefore, W 
quantifies the net quantized circulation of the conden- 
sate around the ring. Thus, breaking of U{1) symmetry 
leads to independent phase selection and results in a net 
superfiuid circulation in the ring [8,]. 

We simulated Eq. ([T]) numerically to study the quench 
dynamics. Growth of particle number A/", total angu- 
lar momentum L and specific angular momentum L/N 
are summarized in Fig. [21 Though the flow involved more 
mass as e decreases, W sticks to a stable value acquired 
right after the symmetry-breaking. Both L and grow 



as long as e decreases, but L/N stabilizes to a steady 
value along with W. This is analogous to the build up of 
persistent current in superfiuid He [1, [2l| . For spatially 
uniform density we have L/N = as observed at the 
end of the quench (last row), when spatial density fluc- 
tuations are mostly ironed out. Final kinetic energy of 
the BEC depends on the final steady value of W (rather 
than on T and final e), and scales with W quadratically. 
W stabilizes (first row) right at the wake of the conden- 
sation, when A^ (second row) is still negligible. Thus, W 
retains phase information from soon after BEC forma- 
tion at the time of symmetry breaking. Once stabilized, 
W is resilient to dissipation and typical ambient thermal 
fluctuations, due to its quantized nature. 

The scaling in Eq. ([6]) is verified by averaging nu- 
merically over > 10^ realizations. The simulation results 
compare favorably with the KZM prediction. Fig. [3fa). 
The square-root dependence of {\W\) on C, Eq.(j6]) is also 
verified there. The values of {\W\) obtained for C = 120 
is twice of that obtained for C = 30 (Fig. [3^). A direct 
comparison between the numerical results and KZM for- 
mula for our model (using [18]) yields, a{W){KZM) = 
5.94, 4.86 and 3.65 versus a{W) (numerical) = 2.33, 1.83 
and 1.35 for tq — 0.01,0.05 and 0.5 respectively. With 
"naive" KZM (choices of phases over ^-sized regions are 
completely independent), mismatch of this order is con- 
sistent with (actually, less than) past experience [2I. [l6|. 

Accurate initial thermalization is crucial in producing 
correct final state. Noise should randomize phase along 
the entire circumference C . This is in effect a diffussive 
process that requires time ~ . This renders accurate 
reproduction of the scaling behavior for large C difficult. 

Adiabatic-impulse transition is at the heart of the scal- 
ing laws predicted by KZM. Here we observe for the first 
time, its manifestation as a scaling-law for the length 
of the impulse reaction time (effective t\ see, e.g., [22]) 
for the condensate density growth. This is an important 
step, since density evolution in a condensation process is 
easier to study experimentally than the associated defect 
dynamics. We compare (Fig. [21 second row) the growth 
of A", with the instantaneous equilibrium value of A" (thin 
blue straight-line segments), obtained from the relation 
e = -~gN/C (valid for smah ~g) 0. Initially, e > 0, and 
A^ is negligible. But after crossing tc, the instantaneous 
equilibrium value of A^ increases linearly with the rate 
N = C/qtq = 0.6 X 10^ tin t = 500rQ, where the ramp 
is stopped and equilibrium value for A" settles to 3 x 10^. 
But for the actual evolution, the length of the period 
from t = = up to the point denoted by A (a sharp 
knee) in respective figures (see also Fig. [Sb), should be 
proportional (and of the order of) t and - in our SGPE 



case - should scale as r, 
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Fig. [3)3 confirms significance of t for density evolution. 
W and L/N also stabilize around this instant A (Fig. [2|). 
Our simulation results confirm this scaling, as shown by 
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FIG. 2: Time evolution of the winding number W, total number of atoms N, total angular momentum L and specific angular momentum 
L/N for random realizations with different tq. For each tq, different quantities for same realizations are plotted in same color. Each 
quench starts from e(— lOOrg) = 100 and continues up to e(t = SOOrg) = —500, where it is held fixed for the rest of the time. The 
black vertical lines indicates tc = (continuous, middle one) and the two on either sides (dashed) indicate ±£. The second row shows 
the adiabatic-impulse transition in BEC growth. While i gives the KZM estimation for the impulse to adiabatic transition boundary, the 
instant denoted by A in the figure marks the point where N shoots up sharply. The thin blue lines indicate the equilibrium values of N 
for the corresponding instantaneous values of e. 



the overlap of N around the knee, when plotted against 
t/t for different tq (Fig. [8)3). After this knee point, N 
catches up rapidly with its equilibrium value. 

Coming to the experiments, we note, unlike other BEC 
relics of symmetry breaking (like solitons), which are dif- 
ficult to resolve due to thermal noise, and decay due to 
dissipation, W bears a very stable and readable signature 
of the underlying phase transition due to its topologi- 
cal stability and integer nature. Statistics of W can be 
presumably studied, e.g, within the experimental setups 
such as [gI, [iqI - [i2| . Experimental study on the growth of 
BEC density in an effective quasi one dimensional ring 
with adjustable cooling rate parameter has been reported 
already by Esslinger group [TTj. They observed the lin- 
ear growth regime, as well as the "knee" feature, that 
provides an experimental counterpart of the effective-t. 
This should allow direct verification of the KZM scal- 
ing law and quantitative determination of the exponent 
for t. Moreover, experimental determination of the ex- 
ponent z may be possible employing the KZM formulae 



z^/2(l + uz) and 1/(1 + i^z) for the measured values of 
and effective-t exponents respectively, and solving 
for z. The scaling law involving N is amenable to more 
accurate experimental determination, since the exponent 
for t scaling is larger than that for the scaling of a{W) 
(1/8 here). Note that for real BEC transition in 3D, 
theoretical prediction is u = 2/3 (close to experimental 
value 0.67 ± 0.13, [23.]) and z = 3/2. 

Now we consider a natural source of experimental 
imperfection that may perturb the scaling behavior, 
namely, inhomogeneity in the potential. An overall non- 
uniformity may have many causes. To be specific, we 
shall think of tilting the torus in the gravitational field. 
This is represented by an additional potential of the form 
¥{1^) = mgsm{a)R[l — cos (i?)] where a represents the 
tilting angle, i9 = x/Ris the angular coordinate denoting 
the position on the ring, R = C/27T and g the gravita- 
tional acceleration. Different stages of condensate forma- 
tion in the tilted ring is shown in Fig. [1] (lower row). Fig- 
ured (a) shows suppression of W as the function of a. In- 
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FIG. 3: (a) The scaling for (|Vl^|> = ^j2Pna{yV) (averaged over 
> 10^ realizations) with tq for C = 30 (open circle) and 120 (filled 
diamonds). Best fits (red line) give exponents 0.1279 ± 0.0034 for 
the C = 30 and 0.1262 ± 0.0075 for C = 120, compared to the 
expected iy/2{l + uz) = 1/8 = 0.125 (Eq. [6]). (|1^|> - VC, i.e., 
ln({|M/^|>C=i2o) - ln({|H^|>c=3o) ~ ln2 (Eq|2]) is also verified (the 
difference shown in the Fig. is ~ 0.71 compared to In 2 0.69). 
The inset shows histograms for distribution of W. (b) A^" vs t/i. 
The "knees", where shoots up (point A), showing perfect overlap 
for different rg's (values marked in the Fig.). 



set of Fig. m^b) shows that the scahng behavior stih per- 
sists for very smah tilting of a = 0.5°, but the exponent 
increases. The main Fig ^p) shows the disappearance of 
the scahng behavior for shghtly bigger tilting {a = 1°). 
Suppression of W could be caused in the tilted ring due 
to the competition between the finite velocity vf of the 
critical front (determined by the critical condition set lo- 
cally by the inhomogeneity [2]) and the sound velocity Vs 
at which the correlation is established between the con- 
densate phase. Formation of phase-inhomogeneity (local 
fiow) is suppressed wherever Vs > vf^ since phase correla- 
tion is maintained throughout the condensation process 
there: the newborn condensate selects the same phase 
as the existing condensate due to local free-energy min- 
imization [H, [25. This implies steeper fall of with 
Tq [2]. But in our case, the above condition is not met 
for a = 0.5°. The same criteria also cannot explain the 
suppression of VF for a = 1° (Fig. IJJd). One possible 
explanation of this suppression may lie in the dissipation 
of kinetic energy prior to the formation of the complete 
BEG ring. The local fiows generated by KZM at the 
bottom of the ring are susceptible to dissipation before 
the topological constraint (that protects the quantized 
circulation) is imposed, i.e., before the BEC is formed 
completely up to the ring top. Such early energy loss 
might leave the condensate with kinetic energy insuffi- 
cient for quantized circulation. We note that tilting a 
BEC ring experimentally and applying that to offset in- 
advertent horizontal misalignment has been achieved re- 
cently [12| , so while we have not developed the theory of 
various imperfections (exemplified here by tilting in the 
gravitational field), control with the accuracies that may 
compensate for (or avoid) such problems seems possible. 

To summarize, we showed that temperature quench 
through the critical point can produce spontaneous circu- 
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FIG. 4: (a) Winding numbers on tilted rings as function of tilting 
angle a for quenching with tq =0.1 (other parameters same as 
Fig. [2|). (b) Suppression of scaling behavior of {|VF|) with tq for 
a = 1° . The inset shows persistence of scaling behavior for smaller 
tilt of a = 0.5° with a larger exponent 0.158 ± 0.011 



lation of BEC in a ring and verified long standing scaling 
predictions of KZM on this. Our demonstration involv- 
ing winding number and condensate density paves a way 
for experimental verification of KZM scaling laws. It also 
provides prescription for experimental determination of 
the critical exponents of the underlying BEC transition. 
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